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Abstract 

We present a study of D = 4 supersymmetric Yang-Mills matrix models with 50(3) 
mass terms based on the Monte Carlo method. 

In the bosonic models we show the existence of an exotic first /second order transition 
from a phase with a well defined background geometry (the fuzzy sphere) to a phase with 
commuting matrices with no geometry in the sense of Connes. At the transition point the 
sphere expands abruptly to infinite size then it evaporates as we increase the temperature 
(the gauge coupling constant). The transition looks first order due to the discontinuity in 
the action whereas it looks second order due to the divergent peak in the specific heat. 

The fuzzy sphere is stable for the supersymmetric models in the sense that the bosonic 
phase transition is turned into a very slow crossover transition. The transition point is 
found to scale to zero with N. We conjecture that the transition from the background 
sphere to the phase of commuting matrices is associated with spontaneous supersymmetry 
breaking. 

The eigenvalues distribution of any of the bosonic matrices in the matrix phase is 
found to be given by a non-polynomial law obtained from the fact that the joint proba- 
bility distribution of the four matrices is uniform inside a solid ball with radius R. The 
eigenvalues of the gauge field on the background geometry are also found to be distributed 
according to this non-polynomial law. 

1 Introduction 

Reduced Yang-Mills theories play a central role in the nonperturbative definitions 
of M-theory and superstrings. The BFSS (Banks-Fischler-Shenker-Susskind) conjecture 
PQ relates discrete light-cone quantization (DLCQ) of M— theory to the theory of iV 
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coincident DO branes which at low energy (small velocities and/or string coupling) is the 
reduction to + 1 dimension of the 10 dimensional U(N) supersymmetric Yang-Mills 
gauge theory [2j. The BFSS model is therefore a Yang-Mills quantum mechanics which 
is supposed to be the UV completion of 11 dimensional supergravity. As it turns out the 
BFSS action is nothing else but the regularization of the supermembrane action in the 
light cone gauge [3]. 

The BMN model [4J is a generalization of the BFSS model to curved backgrounds. It is 
obtained by adding to the BFSS action a one-parameter mass deformation corresponding 
to the maximally supersymmetric pp-wave background of 11 dimensional supergravity. 
See for example [5HZ]- We note, in passing, that all maximally supersymmetric pp-wave 
geometries can arise as Penrose limits of AdS p x S q spaces [8]. 

The IKKT model [9] is, on the other hand, a Yang-Mills matrix model obtained by 
dimensionally reducing 10 dimensional U(N) supersymmetric Yang-Mills gauge theory to 
+ dimensions. The IKKT model is postulated to provide a constructive definition of 
type II B superstring theory and for this reason it is also called type IIB matrix model. 
Supersymmetric analogue of the IKKT model also exists in dimensions d = 3, 4 and 6 
while the partition functions converge only in dimensions d = 4, 6 |48y49j . 

The IKKT Yang-Mills matrix models can be thought of as continuum Eguchi-Kawai 
reduced models as opposed to the usual lattice Eguchi-Kawai reduced model formulated 
in [11]. We point out here the similarity between the conjecture that the lattice Eguchi- 
Kawai reduced model allows us to recover the full gauge theory in the large N theory and 
the conjecture that the IKKT matrix model allows us to recover type II B superstring. 

The relation between the BFSS Yang-Mills quantum mechanics and the IKKT Yang- 
Mills matrix model is discussed at length in the seminal paper [10] where it is also shown 
that toroidal compactification of the D-instanton action (the bosonic part of the IKKT 
action) yields, in a very natural way, a noncommutative Yang-Mills theory on a dual 
noncommutative torus [22]. From the other hand, we can easily check that the ground 
state of the D-instanton action is given by commuting matrices which can be diagonalized 
simultaneously with the eigenvalues giving the coordinates of the D-branes. Thus at tree- 
level an ordinary spacetime emerges from the bosonic truncation of the IKKT action while 
higher order quantum corrections will define a noncommutative spacetime. 

In summary, Yang-Mills matrix models which provide a constructive definition of string 
theories will naturally lead to emergent geometry [18] and non-commutative gauge theory 
[19 , 20j . Furthermore, non- commutative geometry [121,23] and their non- commutative field 
theories [24, 25J play an essential role in the non-perturbative dynamics of superstrings 
and M-theory. Thus the connections between non-commutative field theories, emergent 
geometry and matrix models from one side and string theory from the other side run deep. 

It seems therefore natural that Yang-Mills matrix models provide a non-perturbative 
framework for emergent spacetime geometry and non- commutative gauge theories. Since 
non-commutativity is the only extension which preserves maximal supersymmetry, we 
also hope that Yang-Mills matrix models will provide a regularization which preserves 
supersymmetry [26J. 

In this article we will explore in particular the possibility of using IKKT Yang-Mills 
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matrix models in dimensions 4 and 3 to provide a non-perturbative definition of emergent 
spacetime geometry, non-commutative gauge theory and supersymmetry in two dimen- 
sions. From our perspective in this article, the phase of commuting matrices has no 
geometry in the sense of Connes and thus we need to modify the models so that a geom- 
etry with a well defined spectral triple can also emerge alongside the phase of commuting 
matrices. 

There are two solutions to this problem. The first solution is given by adding mass 
deformations which preserve supersymmetry to the flat IKKT Yang-Mills matrix models 
|12j or alternatively by an Eguchi-Kawai reduction of the mass deformed BFSS Yang-Mills 
quantum mechanics constructed in [13, 60-62] . The second solution, which we have also 
considered in this article, is given by deforming the flat Yang-Mills matrix model in D = 4 
using the powerful formalism of cohomological Yang-Mills theory [30- 32 , 63] . 

These mass deformed or cohomologically deformed IKKT Yang-Mills matrix models 
are the analogue of the BMN model and they typically include a Myers term [14] and 
thus they will sustain the geometry of the fuzzy sphere [27,28] as a ground state which at 
large N will approach the geometry of the ordinary sphere, the ordinary plane or the non- 
commutative plane depending on the scaling limit. Thus a non-perturbative formulation 
of non-commutative gauge theory in two dimensions can be captured rigorously within 
these models [ToTlT?] . See also [64"]l65]. 

This can in principle be generalized to other fuzzy spaces [29] and higher dimensional 
non-commutative gauge theories by considering appropriate mass deformations of the flat 
IKKT Yang-Mills matrix models. 

The problem or virtue of this construction, depending on the perspective, is that in 
these Yang-Mills matrix models the geometry of the fuzzy sphere collapses under quantum 
fluctuations into the phase of commuting matrices. Equivalently, it is seen that the geom- 
etry of the fuzzy sphere emerges from the dynamics of a random matrix theory [521154] . 
Supersymmetry is naturally expected to stabilize the spacetime geometry, and in fact 
the non stability of the non-super symmetric vacuum should have come as no surprise to 
us [33] . 

We should mention here the approach of [34] in which a noncommutative Yang-Mills 
gauge theory on the fuzzy sphere emerges also from the dynamics of a random matrix 
theory. The fuzzy sphere is stable in the sense that the transition to commuting matrices 
is pushed towards infinite gauge coupling at large N [53] . This was achieved by consider- 
ing a very special non-supersymmetric mass deformation which is quartic in the bosonic 
matrices. This construction was extended to a noncommutative gauge theory on the fuzzy 
sphere based on co-adjoint orbits in [35] . 

Let us also note here that the instability and the phase transition discussed here were 
also observed on the non-commutative torus in [39H4T] [69 ] [70] where the twisted Eguchi- 
Kawai model was employed as a non-perturbative regularization of non-commutative 
Yang-Mills gauge theory [36H38] , 

In this article we will then study using the Monte Carlo method the mass deformed 
Yang- Mills matrix model in D = 4 as well as a particular truncation to D = 3. We 
will also derive and study a one-parameter cohomological deformation of the Yang-Mills 
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matrix model which coincides with the mass deformed model in D = 4 when the parameter 
is tuned appropriately. We will show that the first/second order phase transition from 
the fuzzy sphere to the phase of commuting matrices observed in the bosonic models is 
converted in the supersymmetric models into a very slow crossover transition with an 
arbitrary small transition point in the large N limit. We will determine the eigenvalues 
distributions for both D = 4 and D = 3 throughout the phase diagram. 

This article is organized as follows. In section 2 we will derive the mass deformed 
Yang- Mills quantum mechanics from the requirement of supersymmetry and then reduce 
it further to obtain Yang-Mills matrix model in D = 4 dimensions. In section 3 we will 
derive a one-parameter family of cohomologically deformed models and then show that 
the mass deformed model constructed in section 2 can be obtained for a particular value 
of the parameter. In section 4 we report our first Monte Carlo results for the model D = 4 
including the eigenvalues distributions and also comment on the D = 3 model obtained 
by simply setting the fourth matrix to 0. We conclude in section 5 with a comprehensive 
summary of the results and discuss future directions. 

2 Mass Deformation of D = 4 Super Yang-Mills 
Matrix Model 

2.1 Deformed Yang-Mills Quantum Mechanics in 4D 

The M = 1 supersymmetric Yang-Mills theory reduced to one dimension is given by 
the supersymmetric Yang-Mills quantum mechanics (with Dq = do — i[Xq, .]) 

A) = ^TrQpoY,) 2 + l[X t ,X,] 2 - ^ 7 °A)V> + \H\X % M + V) . (2.1) 

The corresponding supersymmetric transformations are 

S X = ej ip 
5 Xi = ejiip 

5 F = -ie757o£>o^ + £757; , Y>] . (2.2) 

Let \i be a constant mass parameter. A mass deformation of the Lagrangian density Cq 
takes the form 

£ M = Co + 4-Ci + + - (2.3) 

The Lagrangian density Co has mass dimension 4. The corrections C\ and C2 must have 
mass dimension 3 and 2 respectively. We recall that the Bosonic matrices Xq and X a have 
mass dimension 1 whereas the Fermionic matrices ipi have mass dimension |. A typical 
term in the Lagrangian densities C\ and C2 will contain rif Fermion matrices, Boson 
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matrices and n t covariant time derivatives. Clearly for C\ we must have ^nj + nb + nt = 3. 
There are only three solutions (n/, n b , nt) = (2, 0, 0), (0, 3, 0), (0, 2, 1). For £2 we must have 
|nj + nb + nt = 2 and we have only one solution (nf,nb,nt) = (0,2,0). Thus the most 
general forms of C\ and £2 are 

A = Tr UMi, + \-S ahc X a X h X c + J ab X a D X b ) . (2.4) 



C2 = Tr(^-^S ab X a X b \ (2.5) 

Clearly for £3 we must have |nj + n b + nt = 1 which can not be satisfied. Thus the 
correction £3 and all other higher order corrections vanish identically. 

We will follow the method of [13] to determine the exact form of the mass deformation. 
We find the fermionic mass term 

£^ = Tri>(iaU + ^;7V,7 j ] + c 7 V7 3 )v<- (2-6) 

The numerical coefficients a, Hij and c will be constrained further under the requirement 
of supersymmetry invariance. 

Next we consider the bosonic terms. By rotational invariance we can choose J ab = 0, 
S a b = v $ab and S a bc = Giee a b c for some numerical coefficients v and e. 

The mass deformed super symmetric transformations will be taken such that on bosonic 
fields they will coincide with the non deformed supersymmetric transformations so that the 
Fierz identity can still be used. The mass deformed supersymmetric transformations on 
fermionic fields will be different from the non deformed supersymmetric transformations 
with a time dependent parameter e = e(t) which satisfies doe = ulle. We will suppose the 
supersymmetric transformations 

5fi,Xi = 5oXi 

W = <W + Z-tAe. (2.7) 



By requiring that the Lagrangian density (|2.3p is invariant under these transformations 
we can determine precisely the form of the mass deformed Lagrangian density and the 
mass deformed supersymmetry transformations. 

A long calculation will yield = 0, c = — 3e/4, v = e 2 - 16a 2 /9, LT = -(2c + 
e)7°7 1 7 2 7 3 , A = — (4ia/3 + e7 1 7 2 7 3 )7*Xj. The mass deformed Lagrangian density and 
mass deformed supersymmetry transformations are given respectively by (with iau = //i/4 
and — 3e^/4 = //2/4) 

£ M = £ + -^ 2 -Tri;(u 1 + / U27 1 7 2 7 3 )^ - ieijk-^TrXiXjXk - ^2(^1 + u l) TrX i. 

(2.8) 
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S^X = e-fQip 

<W = ( - ^[/.fl^o^ + J[f ,7^1^,^] - - ^TS'T 3 )^^)^ (2-9) 

e = e ( t ) = e |(Mi7°-M27V7 2 7 3 )^ ( 2.1 ) 

We verify that {5^X^) + = SuXn and hence the Hermitian matrices Xu remains Hermitian 
under super symmetry. The corresponding supersymmetric algebra is su(2|l) [13J. 

2.2 Truncation to Dimension 

We consider now the Lagrangian density (action) given by 

(2.11) 



A) = -Tr[ -[X„X U ][X^X»} + -^[X^] ) . (2.12) 



^Tr(h,X^X v ][X^X"\ + ^ 

In above we have allowed for the possibility that mass deformations corresponding to the 
reduction to zero and one dimensions can be different by including different coefficients a, 
b and c in front of the fermionic mass term, the Myers term and the bosonic mass term 
respectively. However we will keep the mass deformed supersymmetric transformations 
unchanged. 

After some algebra we find that we must have a = 2/3, c = 1, b = 1 and fj,i = 0. The 
model of interest is therefore 

^T./^y y iivii yi/i i Kt„,h{y „/,i i t i 93/, ^^^2 



ieijk^-TrX^Xk). (2.13) 



Since ip and e are Major ana spinors we can rewrite them as 



*- rv T - rt +)T ■ 



We compute with Xq = 1X4 the action 

^Tr Q[X 4 , X,] 2 + ~ ([XuXj] - i^fe ijk X k j + 9+ [i[X^ ..] + a^X, ..] + ^ 



^iVv. V12 , i/rv. vl .-02 V V , fl +/,ry, 1 , ..r Y . 1 , ^2 



(2.15) 



The supersymmetric transformations are 
5„X = i(w + - 



= i(9 + aiU} - u + o l 



8^9 = (- iai[Xo,Xi] - ^ei jk ak[Xi,Xj} + ^^iXAu. (2.16) 



3 Cohomological Approach 



3.1 Cohomologically Deformed Super symmetry 

The reduction to zero dimension of the M = 1 supersymmetric Yang- Mills theory in 
four dimensions is given by (in Euclidean signature) 

S = --Tr[X^X u f -Trd+(i[X A ,.]+a a [X a ,..] ) jd + 2TrB 2 . (3.1) 

The supersymmetric transformations become 

SX" = icda^e - iOu^uj 

66 = ia' lu [X^,X u }uj -2Buj 

56 = -iOa^iX^, X u ] + 2Buj 

8B = hjjo»[X»,e] + ^[X»,e}a»u. (3.2) 

Let us note that since we are in Euclidean signature the transformation law of X4 is 
antihermitian rather than hermitian. 

By using a contour shifting argument for the Gaussian integral over B we can rewrite 
the auxiliary field B as 

B = H+±[X 1 ,X 2 \. (3.3) 

We will also introduce 

Oi = m + iVi , 62 = Xi + iX2- (3.4) 
cf) = ±{X 3 + iX A ) , $ = -^(X 3 - iX A ). (3.5) 

We compute 
S 

2 '-'cohom 

= Tt(h 2 + H[X 1 ,X 2 ] + [X h <P] [X u 4>] + [</>, 4>f - m [<j>, m ] - xS, Xi] - Vie tJ \Xi, Xj] + %[x», . 

(3.6) 

We have four independent real supersymmetries generated by the four independent grass- 
mannian parameters £j, pi defined by the equations uj\ = £2 + i£i and UJ2 = pi + ip2- 
We look at the supercharge Q\r associated with £2- We define the exterior derivative d 
on bosons by dB = i[Qm,B] and on fermions by dF = i{QiR,F}. The corresponding 
supersymmetric transformations are precisely given by 5B = 2^dB and 5F = 2&dF 
where 



dXi = Xi 
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(3.7) 



d(p = , d(p = —r]2- 



(3.8) 



dH = [4>M- (3-9) 

d Vl = H , d V2 = $,</>]. (3.10) 

d Xi = [<f>,Xi\. (3.11) 

From these transformation laws we can immediately deduce that for any operator O we 
must have 

d 2 = [(/), 0\. (3.12) 

Thus d? is a gauge transformation generated by 4> and as a consequence it is nilpotent on 
gauge invariant quantities such as the action. 
Next we compute 

dTrQ = S cohom . (3.13) 

Q = -Xi[X i ,4>}+r h [X 1 ,X 2 }+r n H- m [ ( f>,4,}. (3.14) 

Thus we have 

d 2 TrQ = dS cohom = 0. (3.15) 

We consider now the deformed action and deformed exterior derivative given by 

Sdef = 5'cohom + S = S'cohom + ^iSi + e 2 S2 + ■■■ (3.16) 

d def = d + eT. (3.17) 

Super symmetric invar iance requires 

ddcfSdcf = 0. (3.18) 

The fact that d 2 is equal on gauge invariant quantities, i.e. d 2 S co \ lom = d 2 Si = leads 
to cPSd c f = 0. We have the identity 

d 2 Ac[ S dc[ = 0. (3.19) 

This is equivalent to 

{d, T}S cohom + eT 2 S cohom + e{d, T}S + eT 2 S = 0. (3.20) 
Thus we must have among other things 

{d,T}S cohom = {d,T}S = 0. (3.21) 



In other words {d, T} generates one of the continuous bosonic symmetries of the action 
Sdef which are gauge transformations and the remaining rotations given by the SO (2) 
subgroup of 50(4). Following [32J we choose {d, T} to be the rotation U defined by 

U : X a — > ie ab X b , Xa — ► ie a bXb- (3.22) 

We have then 

{d,T} = U. (3.23) 

The symmetry T must also satisfy 

T 2 = 0. (3.24) 

By following the method of |44| we can determine precisely the form of the correction T 
from the two requirements T 2 = and {d,T} = U and also from the assumption that T is 
linear in the fields. A straightforward calculation shows that there are two solutions but 
we will only consider here the one which generates mass terms for all the bosonic fields. 
This is given explicitly by 



TXi = , Txi = itijXj , Tcj) = 

TH = i^T} 2 , Trj 2 = , Tr) X = -%\<f> + i 7 , Tcj) = 0. (3.25) 
The cohomologically deformed super symmetric transformations are therefore given by 

d dci Xi = Xi- (3.26) 

d def (f) = , ddef^ = ~m- (3.27) 

d de{ H = [4>,r]i] + ie^r]2- (3.28) 

d&eWi = H + e(-i\<j) + i-/4>) , d dci r] 2 = [<j>, (/>]. (3.29) 

ddeiXi = [(pjXil+ieeijXj. (3.30) 

3.2 Cohomologically Deformed Action 

Next we need to solve the condition (|3.18p . The deformed action is a trace over some 
polynomial P. In the non-deformed case we have S = dQ where Q is a [/—invariant 
expression given by A3. 14j) . We assume that the deformed action S de f = TrP is also 
[/—invariant. By using the theorem of Austing [33] we can conclude that the general 
solution of the condition (|3.18|) . or equivalently of the equation d de fTrP = 0, is 



Sdef = ddefTVQdef + TrR 3 (cj)). 
9 



(3.31) 



For SU(N) gauge group this result holds as long as the degree of P is less than 2N/3. 
Clearly when the deformation is sent to zero cidef — > d, Qdef — > Q and R — > 0. Thus 
we take 



Qdef = Q-iR , R = kiRi + k 2 R 2 . (3.32) 
We choose R\ and R 2 to be the [/—invariant quantities given by 

Ri = ^abXaX b , R 2 = -m<t>- (3-33) 
We choose Rz{<t>) to be the [/—invariant quantity given by 

Rs(4>) = -pV- (3-34) 

In order to remove the deformation we must take e — > so that dd e f — > d and p — > 
so that 5d c f — > dTrQdtf and K{ — > so that Qdef — > Q- 
We compute 

S def = dTrQ - idTrR + eTTrQ - ieTTrR + TrR 3 (<p). (3.35) 

The first term Scohom = dTrQ is the original action. We will choose the parameters so 
that the total action enjoys SO(3) covariance with a Myers (Chern-Simons) term and 
mass terms for all the bosonic and fermionic matrices. We find after some algebra 

7 = , A = 4 - ^ , p 2 = 2e(2e - Kl ). (3.36) 

Thus the total action becomes 

Sdef = •S'cohom + AS^o^m. (3.37) 

AScohom = in 1 Tr{xiX2-r]ir]2) + -eniTrXl- -(k 2 - €j)(k 2 - ej + e\)TrX% 

l o 

- ^(eA + e 7 + 2Ki - K 2 + 4e)TrX 3 [X 1 ,X 2 ] 

= ^Tr6 + e + ^e^TrXl + ± Kl (2e - K X )TrX\ - % -{Ae + Kl )e abc TrX a X b X c . 

(3.38) 

We introduce now 

Qt 

- (4e + ki) = a , Kl = -- + 4C a. (3.39) 



Also we perform the scaling 



A^pAOt^.^/A-L^, ,3.40) 
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and 



2(2iV)4a. 

We get then the one-parameter family of actions given by (we set B = 0) 



(3.41) 



S, 



def 



~Tr[X^, X v \ 2 + N^e abc TrX a X b X c + ^-(1 + 6Co)(l - l2( )TrX 2 a 



1 



+ 4iVa 2 Co(l - l2C )TrXi - —Tr6 + U[X A , ..] + a a [X a , ..} + — (1 - 12^ 



2a 



For stability the parameter Co must be in the range 

< Co < 1/12. 



(3.42) 



(3.43) 



This action for Co = is precisely the mass deformed action derived in section 2. The 
value Co = 1/12 will also be of interest to us in this article. This one-parameter family of 
actions preserves only half of the J\f = 1 supersymmetry in the sense that we can construct 
only two mass deformed supercharges 



4 Simulation Results for D = 4 Yang-Mills Matrix 
Models 

4.1 Models, Supersymmetry and Fuzzy Sphere 

We are interested in the cohomologically deformed Yang- Mills matrix models 
1, 



NTr 



2ia 
Y 



^ [X/j, , X v \ 2 H — — e a b c X a XbX c 



+ Nf3TrX 2 + Nf3 4 TrXl 



2-Tr9+(i[X A , ..] + a a [X a , ..] + c) 9. 



The range of the parameters is 

(3 = |( Q + 6 £ )(a - 12&) , A = 4^ (a - 12ft) , C = |(a - 12ft) 



(4.1) 



(4.2) 



o<c < 



12' 



(4.3) 



This action preserves two supercharges compared to the four supercharges of the original 
non deformed Yang-Mills matrix model [H]. We will be mainly interested in the "mini- 
mally" deformed Yang- Mills matrix model corresponding to the value Co = a/12 for which 
we have 



S 



NTr 



^ [X/j, , X v } 2 -\ — CabcXaXfrXc 



1 

JVof 



Tr9+ (i[X 4 ,..]+a a [X a ,.; 



(4.4) 
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The "maximally" deformed Yang- Mills matrix model corresponding to the value Co = 
coincides precisely with the mass-deformed model in D = 4 and as such it has a full 
M = 1 mass deformed supersymmetry besides the half M = 1 cohomologically deformed 
supersymmetry. From this perspective this case is far more important than the previous 
one. However there is the issue of the convergence of the partition function which we will 
discuss shortly. In any case the "maximally" deformed Yang-Mills matrix model is given 
by the action 



-^[X[i,X u ] 2 -\ — € abc X a X b X c 



S = NTr 

- ^-Tr9 + {i[X 4 ,..]+a a [X a ,..] + |a )(). 



+ N 2 ^TrX 2 a 



The above two actions can also be rewritten as 
SSUSY = NTr 



- 4 P^M' Xv} 2 H ^-CabcXaXbXr 



+ N~Ba 2 TrX 2 a 



- j^Tre+ (i[X A , ..} + a a [X a , ..} + a^j 9. 



(4.5) 



(4.6) 



$ = , £ = cohomologically deformed. 

2 - 2 
B = — , £ = - mass deformed. 
9 3 



(4.7) 



We remark that the bosonic part of the mass-deformed Yang-Mills matrix action can be 
rewritten as a complete square, viz 



S B = NTr 



— -^[Xfj,, X u ] 2 H — —e a b c X a XbX c 



2a 2 9 
+ N—TrX 2 



= NTr(^-[X^X v \ + ^^ uX X x 



(4.8) 



Clearly e^ v \ = if any of the indices takes the value 4. Generically the bosonic 

action of interest is given by 



S B = NTr 



■^[Xf_i, X^] 2 + —^-e abc X a X b X c 



+ N8a 2 TrX 2 . 



(4.9) 



Here we allow 8 to take on any value. The variation of the bosonic action for generic 
values of B reads 

SS B = NTrJ 4 5X 4 + NTrJ b 5X b 

J 4 = [X a , [X a , X 4 ]] , J b = 2Ba 2 X b + i[F ab , X a ] + [Y 4 , [Y 4 , X b }} , 

F ab = i[X a ,X b ] + ae abc X c . (4.10) 

Thus extrema of the model are given by 1) reducible representations J a of SU(2), i.e 
X a = J a an d Y 4 = and 2) commuting matrices, i.e X^ belong to the Cartan sub- 
algebra of SU (N) . The identity matrix corresponds to an uncoupled mode and thus we 
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have SU(N) instead of U(N). Global minima are given by irreducible representations of 
SU(2) of dimensions N and 0. Indeed we find that the configurations X a = <pL a , X4 = 
solve the equations of motion with <\> satisfying the cubic equation (f){4> 2 — a<p + /3a 2 ) = 0. 
We get the solutions 

1 ± J 1-4/3 

00 = 0, 0± = a V — . (4.11) 

We can immediately see that we must have /3<l/4 which does indeed hold for the values 
of interest /3 = and f3 = 2/9. However the action at 4>± is given by 

S fl [fe] = ^c>(/J-l T iVT^). (4.12) 

We can verify that Sb[4>-\ is always positive while Sb[<P+] is negative for the values of 
/3 such that (3<2/9. Furthermore we note that S[<fio\ = 0. In other words for /3<2/9 
the global minimum of the model is the irreducible representation of SU(2) of maximum 
dimension N whereas for /3 > 2/9 the global minimum of the model is the irreducible 
representation of SU(2) of minimum dimension 0. 

At ft = 2/9 we get (f> + = 2a/3 and S[(f> + ] = 0. Thus the configuration X a = ^fL a 
becomes degenerate with the configuration X a = 0. However there is an entire SU(N) 
manifold of configurations X a = 2 ^UL a U + which are equivalent to the fuzzy sphere 
configuration. In other words the fuzzy sphere configuration is still favored although now 
due to entropy. Thus there is a first order transition at ft = 2/9 when the classical ground 
state switches from X a = |L a to X a = as we increase /3 through the critical value 
(3 = 2/9. The two values of interest /3 = and (3 = 2/9 both lie in the regime where the 
fuzzy sphere is the stable classical ground state. 

This discussion holds also for the full bosonic model in which we include a mass term 
for the matrix X4. Quantum correction, i.e. the inclusion of fermions, are expected to 
alter significantly this picture. 

Towards the commutative limit we rewrite the action into the form (with X^ = aD^, 
a = aVN = 1/g, a 4 = 1/g 2 and F ab = i[D a ,D b ] + e abc 4>D c ) 



S = 



^TrFl - ^Tr [e abc F ab D c + + («/>(</> - 1) + ft) TrD\ 



2^N Tr[Da ' D4]2 ~ ^ Tr ^ + (^,..] + * a [D a ,..} +i)0. (4.13) 
The 3rd terms actually cancel for all values of j3. Thus 

S = ^TrF! b -^Tr[e abc F ab D c + 4>Dl]-^Tr[D a ,D 4 ] 2 

- ^Tre + {i[D 4 ,..]+a a [D a ,..]+i)e. (4.14) 

The commutative limit iV — > 00 is then obvious. We write D a = <p(L a + A a ) and we 
obtain 
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(4.15) 

4.2 Path Integral, Convergence and Observables 

In the quantum theory we will integrate over N x N bosonic matrices X^ and N x N 
fermionic matrices 6 + and 9 a . The trace parts of X^, 6>+ and 6 a will be removed since they 
correspond to free degrees of freedom. The partition function of the model is therefore 
given by 

Z = J VX^ V6 V6 + SftrX^) 5(Tr6+) S(Tr6 a ) e " SsuSY 

= J VX^b (TrX^) detV e - SB . (4.16) 



detP = J d9d9 + 5(Tr9 a )5(Tr0+)e^ Trd+V9 . (4.17) 

The integration over the fermions yielded the determinant of the 2(N 2 — 1) x 2(iV 2 — 1) 
dimensional matrix D = i\X±, ..] + a a [X a , ..] + £<x This determinant is positive definite 
since every eigenvalue A of 2? is doubly degenerate |42| . The reason lies in the fact that 
the Dirac operator V = iX^ — iXf- + o a X a — u a X^ + £a is symmetric under the exchange 
of left and right operators, viz under X a o —X^ . A much cleaner proof goes as follows. 
Let \J' be an eigenstate of T> with eigenvalue A, in other words 

i[X^% + o a [X a ,V]+ianl> = XSt. (4.18) 

Taking the hermitian conjugate of this equation we get 

iVU, (*+) T ] - o T a [X a , (*+f ] + eo(*+) r = A(*+f . (4.19) 

In above (\E ,+ ) T is a column vector with components given by Multiplying the above 
equation by c"2 and defining the spinor ^ = o"2(^ /+ ) T we arrive at the equation 

z[X 4 ,tf] + a a [X , = A*. (4.20) 

We have also used the identity a a = — c^cxjo^- We conclude that ^> is also an eigenstate 
of T> with the same eigenvalue A. The spinors \P and \£ are charge conjugate to each 
other. In above we have assumed that A is real which follows from the fact that the Dirac 
operator D = i[X&, ..] + a a [X a , ..] + £a is hermitian. This establishes that the determinant 
det V is positive definite for any configuration X^ and hence the model can be accessed 
directly by Monte Carlo simulation. 

Let us also add that the Dirac operator T> admits an approximate chirality operator 
and hence there is an approximate chiral symmetry in this model beside exact rotational 
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invariance, exact gauge invariance and exact charge conjugation. The existence of chiral 
symmetry (though approximate) means that there should exist more structure in the low 
energy fermionic spectrum. 

The partition function Z is also invariant under the translation — > + eX^ 
where e is a small parameter. Under this coordinate transformation the measure dX^ 
changes to (1 + 4(iV 2 — l)e)dX fl . The bosonic action Sb = S4 + S3 + S% changes to 
Sb + e(454 + 35*3 + 252) under X^ — > X^ + eX^, The determinant, on the other hand, 
changes under X^ — > X^ + eX^ as 

detP — > (l + e) 2(Ar2 - 1) det'^l - -^ 7 ea(l + 7))det'x'. (4.21) 



The matrices Ai and 7 are given in appendix B. We obtain then 

1 z - ' 1 



detP — > 1 + e 



2 (N 2 - 1) - |«Tr ad - - e«Tr ad ^ 7 



JdetP. (4.22) 



From the invariance of the partition function under the coordinate transformation X^ — > 
X^ + eX^ we derive therefore the Schwinger-Dyson identity 

< YM > < YM > < CS > < RAD > ~ COND 6 
= 6. (4.23) 

This is an exact result. 

The operators YM, YMo and CS are the actions given by 

YM = -^Tr[X a ,X b } 2 , YM = -^-Tr[X , X a } 2 , CS = Z J^e abc Tr[X a , X b ]X c . 

(4.24) 

The action RAD is related to the radius of the sphere. It is given by 

RAD = N(3a 2 TrX 2 . (4.25) 
We will define the radius r of the sphere through the relation 

The total bosonic action is given by S = YM + YM + CS + RAD. 

Convergence of Yang-Mills path inetgrals such as the one given by (|4.16p was studied 
extensively in [H] and in |45H50| . This question is of paramount importance for analytical 
analysis as well as for Monte Carlo simulation. The source of the divergence, if any lies 
in the so-called flat directions, i.e. the set of commuting matrices. 

The path integral (|4. 16[) corresponds to a gauge theory with gauge group SU(N) in 
dimension D = 4. We will also consider SU(N) gauge theory in dimension D = 3. We 
start the discussion with the model a = 0, f3 = and £ = 0. It was found in [45] that the 
bosonic path integral in D = 3 is convergent for > 4 while the bosonic path integral in 
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D = 4 is convergent for N > 3. Since we are interested in large values of N we can safely 
consider the bosonic path integrals in D = 3, 4 to be convergent for all practical purposes. 
On the other hand, it is found in [49], that the supersymmetric path integral in D = 3 
is not convergent while the supersymmetric path integral in D = 4 is convergent for all 
N > 2. 

Tuning the parameters a, f3 does not change this picture. For example it was shown 
in [50] that adding a Myers term, i.e. considering a non-zero value of a, does not change 
the convergence properties of the D— dimensional Yang-Mills matrix path integral. The 
point is that the Chern-Simons (Myers) term is always small compared to the quartic 
Yang-Mills term. The same argument should then lead to the conclusion that adding a 
bosonic mass term, i.e. if we consider a non-zero /3, will not change the above picture. 

Tuning the fermion mass term, i.e. considering a non-zero value of the scalar curvature 
£, will lead to complications. In this case the Pfaffian, or equivalently the determinant, 
will be expanded as a polynomial in the scalar curvature £. The analysis of [19] should 
then be repeated for every term in this expansion. We claim that the supersymmetric 
path integral in D = 4 is not convergent for generic values of £. 

We have extensively checked in Monte Carlo simulation the conjecture that Yang-Mills 
matrix models in dimension D = 4 does not make sense for generic values of £. The major 
observation is that for the fermion determinant for generic values of a = y/N a never 
reaches thermalization |^J. However, we have also observed that for sufficiently small values 
of £ the theory actually makes sense and thus there is some critical value of £, which we 
will not determine in this article, above which the path integral is ill defined. The value 
of interest £ = 2/3 corresponding to the mass deformed matrix model lies in this range 
where the model is actually ill defined. 

Therefore in order to access the mass deformed Yang-Mills matrix model by the Monte 
Carlo method we must regularize the theory in such a way as to make sure that the path 
integral is absolutely convergent. Unfortunately most regularizations will not maintain 
neither the full N = 1 mass deformed supersymmetry nor the half M = 1 cohomologically 
deformed supersymmetry of this model. We adopt here the regularization in which we 
will simply set £ = 0. In other words we make the replacement 



•SsiJSY 



S, 



SUSY 



NTr 



—Tr9 + (i[X^..]+a a [X a ,.. 



+ NPa 2 TrX 2 a 



(4.27) 



In summary, the value f3 = 2/9 corresponds to the mass deformed Yang-Mills matrix model 
with softly broken supersymmetry whereas the value j3 = is precisely the minimally 
deformed model which enjoys half of the N = 1 cohomologically deformed supersymmetry. 



1 This happens typically for small values of a far from the fuzzy sphere region but not very close to a — 0. 
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4.3 Bosonic Theory: Emergent Geometry and Phase Dia- 
gram 

Emergent Geometry We measure the different observables as a function of the cou- 
pling constant a for the two relevant values of f3, i.e. f3 = 0,2/9. The first observable is 
the bosonic Schwinger-Dyson equation given by 



TDF - 1 <YM> <YM > q <CS> < RAD > 4 

We have verified that the bosonic Schwinger-Dyson equation holds in Monte Carlo simu- 
lations to a very satisfactory accuracy. 

The radius which we have defined by the equation < 1/r >=< TrX 2 > /a 2 C2 is 
shown on figures and ([2]). For large values of a the result is consistent with the 
classical prediction 



1 TrX 2 „ 1 + a/1 -4/3 

< - >=< 4^ >= <t> 2 + , 0+ = ^ -• (4.29) 

r a C2 2 

This means in particular that the system is in the ground state configurations 



X A = , X a = acf>L a . (4.30) 

In other words we have a fuzzy spherical geometry given by the commutation relations 

[X 4 , X a ] = 0, [X a ,X b ] = ie abc a<pX c . (4.31) 

We have checked these commutation relations and found them to hold quite well for 
sufficiently large values of a. The coordinates on the sphere are defined by 

n„. = 4^, I> 2 a = ^ 2 . (4.32) 



C 2 a 

We observe that as we decrease a, the radius 1/r jumps abruptly to then starts to 
increase again until it becomes infinite at a = 0. This is the most dramatic effect of the 
so-called sphere-to-matrix transition in which the sphere suddenly expands and evaporates 
at the transition points then it starts shrinking to zero rapidly as we lower the coupling 
further. 

This is the interpretation advocated in [51H53] for a similar phenomena observed in 
the case of three dimensional bosonic models. As far as we know this phenomena was 
observed in Monte Carlo simulation first in [54] and it was found in analytical work on 
perturbative three dimensional bosonic models in [S3] and then in |56| . 

The transitions for the bosonic mass deformed model with /3 = 2/9 and the bosonic 
cohomological model with f3 = are observed to occur at the following estimated values 



d* = 4.9 ±0.1 , (3 = 2/9. 
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(4.33) 



q„ = 2.55 ±0.1 ,/3 = Q. (4.34) 

The fuzzy sphere phase corresponds to the region a > a* whereas the matrix phase 
corresponds to the region a < d*. In other words the sphere becomes more stable as we 
make (3 smaller (see below). 

The order of the sphere-to-matrix transition is very difficult to determine. Since the 
ground state configurations are X4 = and X a = a<pL a , the theoretical analysis based 
on the effective potential of the three dimensional model done in |5ip52j should also hold 
here largely unchanged (see below). As a consequence we will only summarize here the 
main points omitting much technical details. 

The specific heat C v = (< S 2 > - < S > 2 )/(N 2 - 1) shown on figure ([3]) diverges 
from the side of the fuzzy sphere with a critical exponent equal 1/2. It is equal to 3 in this 
phase where 1/2 is due to the 2 dimensional U{1) gauge field on the sphere, 1/2 is due to 
the normal scalar field on the sphere and 1/2 is due to the scalar field X4. This critical 
behavior is typical of a second order transition. In the matrix phase the specific heat is a 
constant right up to the transition point and it is equal 1 where each matrix contributes 
1/4. There is no divergence from this side and the critical exponent is 0. In other words 
the behavior above and below the critical coupling are different, which is quite unusual, 
but still from the specific heat we qualify this transition as second order. 

The expectation values of the Yang-Mills action and the Myers (or Chern-Simons) 
action are shown on and ([5]). The expectation values of the total bosonic action for 
the two cases /3 = 2/9 and f3 = are shown on figure ([6]). Prom these observables we 
observe a discontinuity at the transition point. Thus the sphere-to-matrix transition is 
associated with a latent heat equal to A < S >=< S > m atrix — < S > sp here which is 
typical of a first order phase transition. It is straightforward to estimate the value of 
this latent heat. The latent heat is released by going from the matrix phase to the fuzzy 
sphere phase for f3 = whereas for (3 = 2/9 the latent heat is released by going in the 
other direction from the fuzzy sphere phase to the matrix phase. 

As we will see from the discussion of the eigenvalues distributions the matrices in 
the matrix phase are commuting matrices centered around 0. 

Phase Diagram The last point we would like to address within the context of the 
bosonic theory is the construction of the phase diagram in the plane a — (3. We have 
already measured two points of this phase diagram which correspond to the two values 
f3 = 2/9 and f3 = 0. In order to map the phase boundary between the fuzzy sphere and 
the matrix phase we choose other values of (3 and measure for each one of them the critical 
value of a from the discontinuity in the radius 1/r. 

The effective potential of the 4 dimensional bosonic Yang-Mills matrix model in the 
Feynman-'t Hooft background field gauge in the ground state configurations X4 = and 
X a = acj)L a can be calculated using the method of [55]. We find 

V c fi ~4 

= a 

2c 2 



6 4 6 3 ~6 2 ] 
4 3 P 2 



+ 21og0 2 . (4.35) 



The difference with the three dimensional bosonic Yang-Mills matrix model lies in the 
factor of 2 multiplying the logarithm. The critical line can then be obtained following the 
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method of [51]. We get 



0* = |(1 + \ l ~ ^) , at = (4.36) 

This prediction is in a very reasonable agreement with the Monte Carlo data. 



4.4 Dynamical Fermions: Impact of Supersymmetry 

In this section we will discuss the effect of the fermion determinant. First we note that 
simulations with fermions are much more harder than pure bosonic simulations. The main 
source of complication is the evaluation of the determinant which is highly non trivial. 
Thermalization is very difficult and as a consequence taking the limit of large iV is not 
so easy even with the use of the Hybrid Monte Carlo algorithm. In the bosonic case we 
could go as large as N = 100 with very decent number of statistics although in this article 
we have only reported data with N up to N = 16. In the fermionic case we will report 
data with N up to N = 10. 

The first thing we have checked is the Schwinger-Dyson identity (|4.23[) . The super- 
symmetric Monte Carlo data agrees well with the prediction 6 as opposed to the bosonic 
data which agrees with the prediction 4. Note that 6 = 4 + 2 where 4 is the number of 
bosonic matrices and 2 is the number of fermionic matrices. 

The most important order parameter with direct significance to the underlying ge- 
ometry is the radius 1/r. See again figures ([I]) and (j2j). We observe that the transition 
sphere-to-matrix observed in the bosonic theory disappeared completely. Again it seems 
here that there is no major difference between the two models with f3 = 2/9 and (3 = 0. 
It is clear from the structure of the action that the theory with a = can not support 
the fuzzy sphere geometry and thus a phase with commuting matrices must still exist. 
However, the transition to the phase of commuting matrices starting from the fuzzy sphere 
phase seems to be a crossover transition not the second/first order behavior observer in 
the bosonic model. This seems to be confirmed by the behavior of the specific heat, the 
Yang-Mills and Myers actions and the total action shown on figures (|3J), ([IE]) and © 
respectively. The jump and critical behavior in the specific heat and the discontinuity in 
the various actions disappeared. 

We have to note here that the observable < TrX\ > /N diverges in the supersymmetric 
theory with a = (3 = [47J. For the mass deformed theory we have (3 = 2/9 and thus 
the observable < TrX\ > /N always exists. We observe on the second graph of figure ([I]) 
that < TrXl > /N increases as we decrease a towards which is consistent with the fact 
that it will diverge in the limit a — > 0. Qualitatively the same phenomena is observed 
for j3 = on the second graph of figure ([5} with more erratic behavior as we decrease 
a towards 0. However in this case we can not infer that < TrX\ > /N exists for all a 
since /3 = although it looks that it does from the data. From this perspective the mass 
deformed model is better than the cohomologically deformed model. 

We have not succeeded in determining precisely the value at which the crossover tran- 
sition occurs but it seems that it depends on N in such a way that it is pushed to smaller 
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values of a as we increase N. From figures (pQ) we can read that the crossover transition 
for ft = 2/9 occurs at a = 3.13, 2.63, 2.38 and 2.13 for N = 4, 6, 8 and N = 10 respectively. 
A simple fit yields the result 

*t = ^,P = 2/9. (4.37) 

The conjecture that the crossover transition occurs at arbitrarily small values of a in the 
large ./V limit is one of the main results of this article. In this way the fuzzy sphere is truly 
stable in the super symmetric theory and does not decay. In any case we are certain that 
the fuzzy sphere in the supersymmetric theory is more stable compared to the bosonic 
theory and the crossover transition to the matrix phase is much slower. This conclusion 
is similar to that of 



4.5 Eigenvalues Distributions 

Bosonic Theory A powerful set of order parameters is given by the eigenvalues dis- 
tributions of the two matrices X$ and X4 . The eigenvalues distribution of the matrix X4 
is qualitatively the same for all values of the coupling constants a. However, the eigen- 
values distribution of the matrix X3 suffers a major change as we go across the transition 
point. In the fuzzy sphere region the eigenvalues distribution of X3 is given by an N— cut 
distribution corresponding to the N eigenvalues — (N — l)/2, (N — l)/2 whereas in the 
matrix phase the eigenvalues distribution of X3 is identical to the eigenvalues distribution 
of X4. 

The eigenvalues distribution p<±(x4) of the matrix X4 is always centered around 0. In 
the fuzzy sphere phase ^4(^4) depends on the coupling constant a. In the matrix phase 
below the critical value the eigenvalues distribution p±{x±) does not depend on a and 
coincides with the eigenvalues distribution of the non deformed model with a = 0. In this 
region the eigenvalues distributions of the matrices X3 and X4 are identical. 

Motivated by the work [57l|58] it was conjectured [66] that the joint eigenvalues dis- 
tribution of d matrices X\, X^-X^ with dynamics given by a reduced Yang-Mills ac- 
tion should be uniform inside a solid ball of some radius R. We have already checked 
that this conjecture works in three dimensions [59J. We will check now that this con- 
jecture holds also true in four dimensions. Let p(x\,X2, X3, X4) be the joint eigenvalues 
distribution of the 4 matrices X\, X2, X3 and X4. We assume that p(x\,x-2, £3, X4) is 
uniform inside a four dimensional ball of radius R. The normalization condition gives 
p(xi, X2, X3, X4) = I/V4 = 2/tt 2 R 4l . We want to compute the eigenvalues distribution of a 
single matrix, say X4, which is induced by integrating out the other three matrices. We 
compute 



-R i-R 

d,X4 / dx3 

-R J~R 



R i-R 
dx2 I dx\ 

-R J-R 



4tt 



R 



-5- / dx4{R z 

6 -I-R 



(4.38) 
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We obtain therefore the eigenvalues distribution 

P^a) = ^(R 2 -xI)1 (4.39) 

This is precisely the fit with a measured value of R for (3 = and a = given by 

R = 1.826 ±0.004. (4.40) 

The above eigenvalues distribution works better for the theory with /3 = 2/9 as shown on 
figure ([5]) with a similar measured value of R given for a = 0.25 by 

R 2/9 = 1.815 ±0.008. (4.41) 

We have found that these two measured values are almost the same throughout the matrix 
phase. 

We emphasize that ^4(^4) is the eigenvalues distribution of the matrix X4 not only in 
the matrix phase but also in the fuzzy sphere phase with a value of R which depends on 
the coupling constant a. We also emphasize that Piix^) is the eigenvalues distribution 
of the matrix X3 in the matrix phase for the bosonic theory for both values (3 = and 
P = 2/9. 

Another non trivial check for this important conjecture is the theoretical prediction of 
the radius 

<^TrXl> = ^. (4.42) 

This means that in the matrix phase the order parameter < jjTrX 2 >, which is related 
to the radius, is constant. Indeed, this is what we see on figures ([I]) and ([2]) with a good 
agreement between the Monte Carlo measurement and the theoretical prediction. The 
observed value of < j^TrX 2 > is slightly below the theoretical prediction for ft = 2/9 
throughout. This is not the case for /3 = where the Monte Carlo measurement starts 
slightly below R 2 /2 and then rises above it as we approach the transition to the fuzzy 
sphere. 

This can potentially be a serious difference between the tow cases (3 = 2/9 and (3 = 0. 
The transition to the sphere in the case of (3 = 2/9 is in the form of an abrupt jump but 
in the case of (3 = there is a slow rise in the matrix phase as we increase a before the 
actual jump. 

The main conclusion of these successful measurements is the fact that the matrices 
X^ in the matrix phase are commuting and thus they are diagonalizable with a joint 
eigenvalues distribution which is uniform inside a ball of dimension R. 



Supersymmetric Theory In the supersymmetric case we found it much easier to 
compute eigenvalues distributions with the value (3 = 2/9 and thus we will only consider 
here the mass deformed model. A sample of the eigenvalues distributions of the mass 
deformed model is shown on figures ([7]) and (|8|). 

Again it is observed that the eigenvalues distribution of A3 in the fuzzy sphere phase is 
given by an A— cut distribution corresponding to the A eigenvalues —(A — l)/2, (A — 

21 



l)/2 whereas in the matrix phase the eigenvalues distributions of X3 is given by ^4(^3) 
with a much larger value of R given for a = 0.25 by 



The eigenvalues distribution of is always centered around given by pi{x^) with a value 
of R which depends on the coupling constant a. This eigenvalues distribution coincides 
with the eigenvalues distribution of X% deep inside the matrix phase below a = 0.25. 

The eigenvalues distribution p^ is therefore universal in the sense that it describes the 
behavior of the eigenvalues of the matrix X4 for all values of a and all values of j3 and 
the behavior of the eigenvalues of the matrix X3 in the matrix phase for all values of f3. 
We note here the difference between this eigenvalues distribution p^ and the eigenvalues 
distribution of the supersymmetric model with a = f3 = |47] . In the latter case the 
distribution extends from —00 to +00 and goes as 1/x 3 for large eigenvalues. It is not 
clear to us at this stage how the two distributions relate to each other. 

The point at which the eigenvalues distributions of X3 and X4 coincide may be taken 
as the measure for the crossover transition point and thus for N = 10 this occurs at a 
point between a = 1 and a = 0.25. 



4.6 Remarks: D = 3 Yang-Mills Matrix Models and Scalar 
Fluctuations 



The D = 3 Yang-Mills matrix models we can immediately consider here can be ob- 
tained from the above D = 4 models by simply setting the fourth matrix X^ to zero. This 
is different from the IKKT supersymmetric Yang-Mills matrix model in D = 3 by the 
fact that it involves a determinant instead of a Pffafian and as a consequence the path 
integrals of these theories are convergent. 

The physics of these D = 3 models is identical to the physics of the D = 4 models in 
the sense that there is a first/second order phase transition from a background geometry 
(the fuzzy sphere) to commuting matrices which in the presence of dynamical fermions 
is turned into a slow crossover transition. The most important difference is the natural 
expectation that the eigenvalues distributions of the matrices X a in the matrix phase must 
be distributed according to the formula 



By analogy with the D = 4 formula (|4.39p this distribution can be derived from the 
assumption that the joint eigenvalues distribution of the 3 matrices X\, X2, X3 is uni- 
form inside a three dimensional ball of radius R. In the next section we will also derive 
this distribution for the D = 4 Yang-Mills matrix model with a particular choice of the 
parameters with (3^ 7^ (see (|4.ip ). 

Monte Carlo measurment of the radius gives the value 



# = 2.85^0.009. 



(4.43) 



3 

^3(^3) = -^(R 2 -x 2 3 ). 



(4.44) 



R ~ 2. 



(4.45) 
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As shown on the first graph of the figure ([9]) a sample of the data for the N = 10 three- 
dimensional bosonic Yang-Mills matrix model with /3 = 2/9 and a = 0.5 is shown. Clearly 
it can be fit nicely to (14.44p with R = 2. In performing the fitting in three dimensions 
we have to cut the tails in order to get a sensible answer. As it turns out the same three 
dimensional data can also be fit to the four dimensional prediction (|4.39|) . We also show 
the data for the TV = 10 four-dimensional bosonic Yang-Mills matrix model with /3 = 2/9 
and a = 0.5 for comparison. 

In the Monte Carlo data of the D = 3 Yang-Mills matrix models reported in this 
article it was not possible to resolve the ambiguity between the two fits (|4.39p and ()4.44p . 
However high precision runs performed in |67| seems to indicate that indeed the three 
dimensional prediction (|4.44[) is the correct behavior for the eigenvalues distributions of 
the D = 3 Yang-Mills matrix models. 

The second remark we would like to discuss in this section concerns the dependence on 
N and a of the eigenvalues distributions p4 and p3 given in (|4.39p and (|4.44p respectively. 
As shown on the second graph of the figure ([9]) the distributions p4 and pj, are independent 
of a. Similarly we can show that these distributions are independent of N . 

The third remark concerns the eigenvalues distributions of the normal scalar field in 
the fuzzy sphere phase which is define by 4> = (X% — 2 C2)/(20y / C2) [55]. The behavior 
in both the D = 3 and D = 4 models is the same although we have to note that the 
effective values C2^s in L> = 3 and D = 4 are slightly different. This is quite natural as 
the three dimensional model is more stable in the sense that it has a lower critical value 
a*. The central observation here is that we can nicely fit these eigenvalues distributions 
to the eigenvalues distribution p± given in (|4.39p . 

5 Summary and Future Directions 

In this article we employed the Monte Carlo method to study nonperturbatively Yang- 
Mills matrix models in D = 4 with mass terms. We can summarize the main results, 
findings and conjectures of this work as follows: 

• By imposing the requirement of supersymmetry and 50(3) covariance we have shown 
that there exists a single mass deformed Yang-Mills quantum mechanics in D = 4 
which preserves all four real supersymmetries of the original theory although in a 
deformed form. This is the 4 dimensional analogue of the 10 dimensional BMN 
model. Full reduction yields a unique mass deformed D = 4 Yang-Mills matrix 
model. This latter 4 dimensional model is the analogue of the 10 dimensional IKKT 
model. 

• By using cohomological deformation of supersymmetry we constructed a one-parameter 
(Co) family of cohomologically deformed D = 4 Yang-Mills matrix models which pre- 
serve two supercharges. The mass deformed model is one limit (Co — > 0) of this 
one-parameter family of cohomologically deformed Yang-Mills models. 

• We studied the models with the values /3 = and f3 = 2/9 where (3 is the mass 
parameter of the bosonic matrices X a . The second model is special in the sense that 
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classically the configurations X a ~ L a , X 4 = is degenerate with the configuration 
X a = 0,X A = 0. 

The Monte Carlo simulation of the bosonic D = 4 Yang-Mills matrix model with 
mass terms shows the existence of an exotic first/second order transition from a 
phase with a well defined background geometry given by the famous fuzzy sphere 
to a phase with commuting matrices with no geometry in the sense of Connes. The 
transition looks first order due to the jump in the action whereas it looks second 
order due to the divergent peak in the specific heat. 

The fuzzy sphere is less stable as we increase the mass term of the bosonic matrices 
X a , i.e. as we increase ft. For ft = 2/9 we find the critical value a* = 4.9 whereas 
for ft = we find the critical value ft = 2.55. 

The measured critical line in the plane a — ft agrees well with the theoretical pre- 
diction coming from the effective potential calculation. 

The order parameter of the transition is given by the inverse radius of the sphere 
defined by 1/r = TrX 2 / (a 2 C2). The radius is equal to l/(fi 2 (where <fi is the classical 
configuration) in the fuzzy sphere phase. At the transition point the sphere expands 
abruptly to infinite size. Then as we decrease the inverse temperature (the inverse 
gauge coupling constant) a, the size of the sphere shrinks fast to 0, i.e. the sphere 
evaporates. 

The fermion determinant is positive definite for all gauge configurations in D = 
4. We have conjectured that the path integral is convergent as long as the scalar 
curvature (the mass of the fermionic matrices) is zero. 

We have simulated the two models ft = and ft = 2/9 with dynamical fermions. The 
model with (3 = has two supercharges while the model ft = 2/9 has a softly broken 
supersymmetry since in this case we needed to set by hand the scalar curvature to 
zero in order to regularize the path integral. 

Thus ft = is amongst the very few models (which we known of) with exact super- 
symmetry which can be probed and accessed with the Monte Carlo method. 

The fuzzy sphere is stable for the supersymmetric D = 4 Yang-Mills matrix model 
with mass terms in the sense that the bosonic phase transition is turned into a very 
slow crossover transition. The transition point a is found to scale to zero with N. 
There is no jump in the action nor a peak in the specific heat. 

The fuzzy sphere is stable also in the sense that the radius is equal 1 /(p 2 over a much 
larger region then it starts to decrease slowly as we decrease the inverse temperature 
a until it reaches at a = 0. We claim that the value where the radius starts 
decreasing becomes smaller as we increase N. 

The model at a = can never sustain the geometry of the fuzzy sphere since it is the 
non deformed model so in some sense the transition to commuting matrices always 
occurs and in the limit N — > 00 it will occur at a* — > 0. 

We have spent a lot of time in trying to determine the eigenvalues distributions of the 
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matrices in both the bosonic and supersymmetric theories. A universal behavior 
seems to emerge with many subtleties. These can be summarized as follows: 

— In the fuzzy sphere the matrices X a are given by the SU{2) irreducible rep- 
resentations L a . For example diagonalizing the matrix X3 gives N eigen- 
values between (N — l)/2 and — (N — l)/2 with a step equal 1, viz m = 
(N - l)/2, (N - 3)/2, -(N - 3)/2,-(N - l)/2. 

— In the matrix phase the matrices X^ become commuting. More explicitly the 
eigenvalues distribution of any of the matrices X^ in the matrix phase is given 
by the non-polynomial law 

P4{X) = 3^ {R2 ~ X2)I - (5 - 1} 

This can be obtained from the conjecture that the joint probability distribution 
of the four matrices X^ is uniform inside a solid ball with radius R. 

— In the matrix phase the eigenvalues distribution of any of the X a , say A3, is 
given by the above non-polynomial law with a radius R independent of a and 
N. 

— This is also confirmed by computing the radius in this distribution and compar- 
ing to the Monte Carlo data. 

— A very precise measurement of the transition point can be made by observing 
the point at which the eigenvalues distribution of A3 undergoes the transition 
from the TV-cut distribution to the above non-polynomial law. 

— The eigenvalues distribution of A4 is always given by the above non-polynomial 
law, i.e. for all values of a, with a radius R which depends on a and N. 

— Another signal that the matrix phase is fully reached is when the eigenvalues 
distribution of A4 coincides with that of A3. From this point downward the 
eigenvalues distribution of X4 ceases to depend on a and N. 

— Monte Carlo measurements seems to indicate that R = 1.8 for bosonic models 
and R = 2.8 for supersymmetric models. The distribution becomes wider in the 
supersymmetric case. 

— We have also observed that the eigenvalues of the normal scalar field X% — &i in 
the fuzzy sphere are also distributed according to the above non-polynomial law. 
This led us to the conjecture that the eigenvalues of the gauge field on the back- 
ground geometry are also distributed according to the above non-polynomial 
law. Recall that the normal scalar field is the normal component of the gauge 
field to the background geometry which is the sphere here. 

• In the D = 3 Yang-Mills matrix model with mass terms the eigenvalues distribution 
becomes polyonomial (parabolic) given by 

P3 (x) = ^(R 2 -x 2 ). (5.2) 

It was difficult for us in this article to differentiate with certainty between the two 
distributions and ps in the three dimensional setting. 
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• Finally, we conjecture that the transition from a background geometry to the phase 
of commuting matrices is associated with spontaneous supersymmetry breaking. In- 
deed mass deformed supersymmetry preserves the fuzzy sphere configuration but 
not diagonal matrices. 

Among the future directions that can be considered we will simply mention the fol- 
lowing four points: 

• Higher precision Monte Carlo simulations of the models studied in this article is the 
first obvious direction for future investigation. The most urgent question (in our 
view) is the precise determination of the behavior of the eigenvalues distributions in 
D = 4 and D = 3. An analytical derivation of P3 and especially p± is an outstanding 
problem. 

• Finding matrix models with emergent 4 dimensional background geometry is also an 
outstanding problem. 

• Models for emergent time, and to a lesser extent emergent gravity, and as a conse- 
quence emergent cosmology are very rare. 

• Monte Carlo simulation of supersymmetry based on matrix models seems to be a 
very promising goal. 

• A complete analytical understanding of the emergent geometry transition observed 
in Yang-Mills matrix models with mass terms is also an outstanding problem. 

In [68], we have attempted to compute the above eigenvalues distributions analyti- 
cally. Using localization techniques we were able to find a special set of parameters 
for which the D = 4 Yang- Mills matrix model with mass terms can be reduced to the 
three dimensional Chern-Simons (CS) matrix model. The saddle-point method leads 
then immediately to the eigenvalues distributions p^. We believe that our theoreti- 
cal prediction for the value of R is reasonable compared to the Monte Carlo value. 
We have also made a preliminary comparison between the dependence of R on a in 
the hermitian and antihermitian CS matrix models. The hermitian case seems more 
appropriate for the description of the eigenvalues of A3 whereas the antihermitian 
case may be relevant to the description of the eigenvalues of X4. 
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Figure 1: The radius of the mass deformed model. 
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1/r of the cohomologically deformed model 
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Figure 2: The radius of the cohomologically deformed model. 



32 



The specific heat of the cohomologically deformed model 
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Figure 3: The specific heat of the supersymmetric models. 
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The YM action of the mass deformed model 
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Figure 5: The average of the Yang-Mills and Myers actions of the mass deformed model. 
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The action of the cohomologically deformed model 
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EVs of the N=10 mass deformed model 
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Figure 7: The eigenvalues distributions of the mass deformed model. 
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EVs of the N=10 mass deformed model with a =0.25 
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EVs of the N=10 mass deformed model with a =0.25 
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Figure 8: The eigenvalues distributions of the mass deformed model. 
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EVs of the N=10 bosonic models with p =2/9 and a =0.5 
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Figure 9: The eigenvalues distributions of X3 of the = 3, 4 bosonic models with (3 = 2/9 in the 
matrix phase. 



39 



